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Abstract 

We derive a scaling property from a fundamental nonlinear differential equation whose solution 
is the so-called (/-exponential function. A scaling property has been believed to be given by a 
power function only, but actually more general expression for the scaling property is found to be 
a solution of the above fundamental nonlinear differential equation. In fact, any power function 
is obtained by restricting the domain of the (/-exponential function appropriately. As similarly 
as the correspondence between the exponential function and Shannon entropy, an appropriate 
generalization of Shannon entropy is expected for the scaling property. Although the (/-exponential 
function is often appeared in the optimal distributions of some one-parameter generalized entropies 
such as Renyi entropy only Tsallis entropy is uniquely derived from the algebra of the (/-exponential 
function, whose uniqueness is shown in the two ways in this paper. 
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I. SCALING PROPERTY DERIVED FROM A FUNDAMENTAL NONLINEAR 
DIFFERENTIAL EQUATION 



The exponential function is often appeared in every scientific field. Among many prop- 
erties of the exponential function, the linear differential function dy/dx = y is the most 
important characterization of the exponential function. A slightly nonlinear generalization 
of this linear differential equation is given by 

Tx =yq (gGM) - (1) 

(See the equation (17) at page 5 of [| and the equations (22)-(23) at page 8 of This 
nonlinear differential equation is equivalent to 
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dy = / dx. (2) 

y q 1 



Then we define the so-called q-logarithm ln„x 



x l ~ q — 1 

K x -= x _ (3) 
as a generalization of lnrr. Applying the property: 

Tx lYi « x = b (4) 

to (J2J), we obtain 

In, y = x + C (5) 

where C is any constant [3]. Then we define the so-called q- exponential exp q (x) as the 
inverse function of ln„ follows: 

( [l + {l-q)x]^ if l + (l-q)x> 0, 
exp g (x) := < (6) 
I otherwise. 

Note that the g-logarithm and g-exponential recover the usual logarithm and exponential 
when q — > 1, respectively (See the pages 84-87 of 0] for the detail properties of these gener- 
alized functions \n q x and exp g (x)). Thus, the general solution to the nonlinear differential 
equation (0) becomes 

V = exp„ (x + C)= exp (C) exp X _ Nl _ q ) (7) 



(exp, (C)) 



where C is any constant satisfying 1 + (1 — q) C > 0. Dividing the both sides by exp q (C) 
of the above solution, we obtain 

eX PQ ~, X „,sl-a ) • ( 8 ) 



Under the following scaling: 
we obtain 

yt = exp q (xf) . (10) 

This means that the solution of the nonlinear differential equation (0) obtained above is 
"scale-invariant" under the above scaling Q. Moreover, we can choose any constant C 
satisfying 1 + (1 — q) C > because C is an integration constant of (J2J). 

Note that the above scaling © with respect to both variables x and y can be observed 
only when q ^ 1 and C ^ 0. In fact, when q — 1, i.e., y = exp (x + C) , © reduces to the 
scaling with respect to only y, i.e., x/ = x, and when C = 0, both scalings in © disappears 

Q 

The above scaling property of the nonlinear differential equation (pQ) is very significant 
in the fundamental formulations for every generalization based on We summarize the 
important points in the above fundamental result. 

1. C in the scaling Q is any constant satisfying 1 + (1 — q) C > because C is an 
integration constant of (J2J). This means that the scaling Q is arbitrary for any q and 
C if q and C satisfies 1 + (1 - q) C > (q ^ 1 and C ^ 0). 

2. In general studies of differential equation, C is determined by the initial condition of 
the nonlinear differential equation ((T}. This means, when an observable in a dynamics 
grows according to the nonlinear differential equation (fl|). i/ie initial condition deter- 
mines the scaling of the dynamics. This is applicable to the analysis of the chaotic 
dynamics SQBBBQH. 

3. In general, for a mapping / : X — > Y, f : Xl (c X) — > Y is called a restriction of a 
mapping / to Xf, which is denoted by / \ Xf. Let f q be a g-exponential function (JBJ) 
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from R — > R + . For the restricted domain R/ g defined by 

», :={xGl|(l-(|)i»l}(cl), (11) 
a restriction of a mapping f q to WL/ q is denoted by 

//, := /, r M/, : M/, -> R + - (12) 

Then // ? becomes a power function: 

f/ q (x) = ((l-q)x)^ (13) 

In the above formulations the only case g < 1 is discussed . As shown in the section IV, 
the ^-generalizations along the line of (JTJ) has a symmetry g <->• 2 — g, i.e., 1 — (1 — g) <-> 
1 + (1 — g). Therefore, the above ff q (x) can be replaced by a restriction of a mapping 
/ 2 _ g to R/ 2 _ 9 : 

// 2 „ 9 (x) = ((g-l)x)^^x^i (14) 

in accordance with the symmetry, which implies that the case g > 1 can be discussed. 
In this way, the restriction of the g-exponential function to the domain M.f q (g < 1) 
or M/2-g (g > 1) coincide with a power function, which has been often appeared and 
discussed in science. In general, the restricted domain W g or R/ 2 -g is called "scaling 
domain" and its corresponding range is called "scaling range" . 

4. A power function / : R — ► R is known to be characterized by the following functional 
equation, i.e., there exists a function g : R — > R such that 

/ (fee) = (6) / (x) (15) 

holds for any 6, x G R. The above functional equation uniquely determines a power 
function 

f(x) = f(l)x- a (16) 

for choosing g (6) = See some references such as for the proof. On the other 
hand, a g-exponential function is characterized by the nonlinear differential equation 
(jTj) as similarly as a exponential function. Moreover, the solutions of are scale- 
invariant under the scaling (0) and reduce to power functions when the domain is 
restricted to W q or R/ 2 -q as shown above. In fact, by restricting the domain of (JSJ) to 
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Wtg, the general solution (jHJ) of the nonlinear differential equation (JTJ) reduces to the 
following power function according to (fH?j) . 



V n \ x \ 9 



TTTn^o = _ 77n (17) 



exp 9 (C) y v fexpJC)) 1 9 J ex P(J (C) 

that is, 

which is equivalent to (fTBj) if 
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(18) 



a=^.. (19) 

Therefore, many discussions on "exponential versus power-law", i.e., = y versus 
f (bx) = g (b) f (x)" should be replaced by "exponential versus g-exponential", i.e., 
"ftc = versus ^ = y qn , which is more natural from mathematical point of view. 

As shown in these discussions, the fundamental nonlinear differential equation {1} pro- 
vides us with not only the characterization of the g-exponential function but also the scaling 
property in its solution. 

As similarly as the relation between the exponential function exp (x) and Shannon en- 
tropy, we expect the corresponding information measure to the g-exponential function 
exp g (x). There exist some candidates such as Renyi entropy, Tsallis entropy and so on. 
But the algebra derived from the g-exponential function uniquely determines Tsallis en- 
tropy as the corresponding information measure. In the following sections of this paper, we 
present the two mathematical results to uniquely determine Tsallis entropy by means of the 
already established formulations such as the g-exponential law, the g-multinomial coefficient 
and g-Stirling's formula. 



II. ^-EXPONENTIAL LAW 



The exponential law plays an important role in mathematics, so this law is also expected 
to be generalized based on the q-exponential function exp (x). For this purpose, the new 



multiplication operation <S> q is introduced in 



ljj and Q 



for satisfying the following identi- 



ties: 



ln 9 (x ® q y) = hiq x + ln 5 y, 
exp 9 (x) <g>, exp g (y) = exp q (x + y) 



(20) 
(21) 
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The concrete form of the g-logarithm or g-exponential has been already given in the previous 
section, so that the above requirements as g-exponential law leads us to the definition of ® g 
between two positive numbers. 



Definition 1 For two positive numbers x and y, the q-product ® q is defined by 

j [x l ~ q + y x - q - l] 1 ^ , if x > 0, y > 0, x 1 ^ + y 1 ^ - 1 > 0, 
x ® q y := < (22) 

I 0, otherwise. 
The q-product recovers the usual product such that lim (x ®„ y) = xy. The fundamental 

q-*l 

properties of the g-product ® q are almost the same as the usual product, but 

a(x® q y) ^ ax® q y (a,x,yEM). (23) 

The other properties of the q-product are available in Q and Q- 

In order to see one of the validities of the g-product, we recall the well known expression 
of the exponential function exp (x) given by 

/ x\ n 

exp (x) = lim 1 1 + - I . (24) 
Replacing the power on the right side of (|2~%j) by the n times of the g-product <g) q : 



(25) 



n times 



exp (x) is obtained. In other words, lim (1 + -) 9 coincides with exp (x) . 

/ X \ ®? 

exp (x) = lim 1 + - (26) 

n— too V n/ 

The proof of ()26|) is given in the appendix of This coincidence ()26|) indicates a validity 
of the g-product. In fact, the present results in the following sections reinforce it. 

III. g-MULTINOMIAL COEFFICIENT AND g-STIRLING'S FORMULA 

We briefly review the g-multinomial coefficient and the g-Stirling's formula by means of 
the g-product ® g . As similarly as for the g-product, q-ratio is introduced as follows: 
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Definition 2 For two positive numbers x and y, the inverse operation to the q-product is 
defined by 



x(Z) q y 



\x 



0. 



1-9 _ yi-q + i] i=5 , if x > 0, y > 0, x 1 " 9 - y l ~ q + 1 > 0, 
otherwise 



(27) 



which is called q-ratio in 

5 is also derived from the following satisfactions, similarly as for ® g Q| • 

ln g (x g y) = h\ q x - \n q y, (28) 
exp g (x) q exp ? (y) = exp q (x - y) . (29) 

The g-product and g-ratio are applied to the definition of the g-multinomial coefficient 

0. 

Definition 3 For n = Xli=i n i an d rii & N (i — 1, ■ ■ ■ , k) , the q-multinomial coefficient is 
defined by 

n 



ni ■■ ■ n k 



'■= ( n] -q) 0g [( n l ] -q) ®g • • • ®g ( n k ] -q)} ■ 



(30) 



J q 



From the definition f)30j) . it is clear that 



lim 

9—1 



n 

m ■ ■■ n k 



n 

m ■ ■ ■ n k 



rv. 



ni! • • -n k \ 



(31) 



In addition to the g-multinomial coefficient, the g-Stirling's formula is useful for many 
applications such as our main results. By means of the g-product (J22J), the g-factorial n\ q 
is naturally defined as follows. 

Definition 4 For a natural number neN and q G R + , the q-factorial n\ q is defined by 

l® q ---® q n. (32) 



rv. 



q ■' 



Using the definition of the g-product (|2*2*|) . ln g (n! 9 ) is explicitly expressed by ln g (n\ q ) = 
^~ ,k= \ k _ -.If an approximation of ln g (n\ g ) is not needed, this explicit form should be di- 
rectly used for its computation. However, in order to clarify the correspondence between 
the studies g = 1 and g ^ 1, the approximation of ln g (n\ q ) is useful. In fact, using the 
following g-Stirling's formula, we obtain the unique generalized entropy corresponding to 
the g-exponential function exp ? (x), shown in the following sections. 



Theorem 5 Let n\ q be the q-factorial defined by Iffjj)). The rough q-Stirling's formula 
ln g {n\ q ) is computed as follows: 



n 



liig (n\ q ) 



ln g n 



n 



+ 0(\n q n) if q^2, 
if q = 2. 

The proof of the above formulas ()33j) is given in [lq . 



2-q H 2-q 
n — Inn + O (1) 



(33) 



IV. TSALLIS ENTROPY UNIQUELY DERIVED FROM THE (/-MULTINOMIAL 
COEFFICIENT AND g-STIRLING'S FORMULA 



In this section we show that Tsallis entropy is uniquely and naturally derived from the 
fundamental formulations presented in the previous section. In order to avoid separate 
discussions on the positivity of the argument in (}3T)j) . we consider the q- logarithm of the 
g-multinomial coefficient to be given by 



ln„ 



n 

n x ■■■ n k 



\n q (n\ q ) - ln g (m! 9 ) ln 9 (n k \ q ) 



(34) 



J q 



The unique generalized entropy corresponding to the g-exponential is derived from the q- 
multinomial coefficient using the g-Stirling's formula as follows jlq . 

Theorem 6 When n is sufficiently large, the q-logarithm of the q-multinomial coefficient 
coincides with Tsallis entropy 1136}) as follows: 

n 



ln„ 



m ■ ■ ■ n k 



where S q is Tsallis entropy 



- n S2-J—,'--,—) if 9>0, q + 2 

2 — q V n n J 

-S 1 {n) + 1 ts 1 (n i ) if q = 2 

i=i 



(35) 



i=l 



S q :— - 

q- 1 



(36) 



and Si (n) is given by Si (n) := Inn. 



The proof of this theorem is given in [15]. 



S 



Note that the above relation (|35|1 reveals a surprising symmetry: ()35|) is equivalent to 



hii- 



(1-9) 



n 

n x •■■ n k 



n k 



+ (1 — q) V n n 



(37) 



J l-(l-g) 

for q > and q ^ 2. This expression represents that there exists a symmetry with a factor 
1 — g around g = 1 in the algebra of the g-product. Substitution of some concrete values of 
q into ()35|) or (pITj) helps us understand the meaning of this symmetry. 

Remark that the above correspondence ()35j) and the symmetry ()37j) reveals that the 
g-exponential function (jUJ) derived from Q is consistent with Tsallis entropy only as in- 
formation measure. 



V. THE GENERALIZED SHANNON ADDITIVITY DERIVED FROM THE q- 
MULTINOMIAL COEFFICIENT 

This section shows another way to uniquely determine the generalized entropy. More pre- 
cisely, the identity derived from the g-multinomial coefficient coincides with the generalized 
Shannon additivity which is the most important axiom for Tsallis entropy. 

Consider a partition of a given natural number n into k groups such as n = Y2^=i n i- 
In addition, each natural number rii (i = 1, • • • , k) is divided into m; groups such as rii = 
Y2T=i ni i w here n%j G N. Then, the following identity holds for the g-multinomial coefficient. 



n 




n ll n 12 • • * n \m x n i\ "' n i mi n kl-" n km k 



FIG. 1: partition of a natural number n 



n 



Tin • • • n h 



n 

ni • • • n k 



- 1 <? 



ni 



n k i ' ' ' n krnk 



(38) 
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It is very easy to prove the above relation ([38)1 by taking the g-logarithm of the both sides 
and using 

On the other hand, the above identity (|38|) is reformed to the generalized Shannon addi- 
tivity in the following way. Taking the g-logarithm of the both sides of the above relation 
()38|) . we have 



ln„ 



n 



n\\ ■ ■ ■ nkm k 
From the relation (J33j) . we obtain 







k 


= ln g 


n 






ni ■ ■ ■ n k 





ri; 



Tin ' ' ' Tli mi 



(39) 



c fn n n kmk 

>~>2-q I , " - - , 

V n n 



S- 



2-q 



->-.-)+E(-fV.(- 

n n/ z — ' \ n I \ n- 

i=l 



- 1 q 



ri; 



Then, by means of the following probabilities defined by 



Pij ■ 



^ (< = 1, 
n 

rrii mi 



k, j = 1, - • • ,m k ) , 



3=1 3=1 



V rii 



£ 

3=1 
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(40) 

(41) 
(42) 



the identity pUjl becomes 

k ( 

Sq (Pll, • • • ,PfemJ = ^ (Pi, ,Pk)+ ^P^Sq (—,••• 

i=l v Pi 



Pi 



Pi 



(43) 



The formula (|43|) obtained from the g-multinomial coefficient is exactly same as the gener- 
alized Shannon additivity (See [GSK3] given below) which is the most important axiom for 



Tsallis entropy 



In fact, the generalized Shannon-Khinchin axioms and the uniqueness theorem for the 
nonextensive entropy are already given and rigorously proved in |l7[. The present result 
(|43|) and the already established axiom [GSK3] perfectly coincide with each other. 



Theorem 7 Let A n be defined by the n- dimensional simplex: 



A, 



(Pl,---,Pn) 



Pi > o, ^Pi = i 



i=l 



The following axioms [GSK1]~[GSK4] determine the function S q : A n 

n 

S q (p l ,...,p n ) = ■ 



(?) 



where 4>(q) satisfies properties (i) ~ (iv): 



(44) 



such that 



(45) 
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(i) (g) is continuous and has the same sign as q — l,i.e., 

0(g)(g-l)>O; (46) 

(ii) 

lim0(g) = 0(l) = O, 0(g)^O(g^l); (47) 

(iii) there exists an interval (a, 6) C R + such that a < 1 < b and (g) is different iable on 
the interval 

(a,l)U(l,6); (48) 

and 

(iv) there exists a constant k > such that 

hm^M = 1 (49 ) 
9->i dq k 

[GSK1] continuity: S q is continuous in A n and g € M + , 

[GSK2] maximality: for any g G R + , any n6N and any (pi, ■ • • ,p n ) G A n , 

Sj(Pl, ••• ,Pn) < S q (-,...,-) , (50) 

[GSK3] generalized Shannon additivity: if 

rrii 

Pij > 0, Pi = ^Pij Vz = 1, • • • , n, Vj = 1, • • • , m h (51) 

3=1 

then the following equality holds: 



(52) 



S q {pil,--- ,Pnm k ) = Sqip!,--- ,p n ) + J^P-^g (—,••• , — ) , 

[GSK4] expandability: 

S x (p l9 . . . ,p„,0) = 5i (pi, . . .,p n ) . (53) 

Note that, in order to uniquely determine the Tsallis entropy (|36|) in the above set of the 
axioms, "lim" should be removed from (|35J). that is, = | (i.e., 0(g) =| (g — 1)) should 
be used instead of ()49|) . The general form 0(g) perfectly corresponds to Tsallis' original 
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introduction of the so-called Tsallis entropy in 1988 [16]. See his original characterization 
shown in page 9 of Q for the detail {(f) (g) corresponds to "a" in his notation. His simplest 
choice of "a" coincides with the simplest form of 0(g) i.e., = |.). 

When one of the authors (H.S.) submitted the paper [17| in 2002, nobody presented the 
idea of the g-product. However, as shown above, the identity on the g-multinomial coefficient 
[15I ] which was formulated based on the g-product ^| coincides with one of the axioms ( 
[GSK3] : generalized Shannon additivity) in 3| ■ This means that the whole theory based on 
the g-product is self-consistent. Moreover, other fundamental applications of the g-product, 
such as law of error Il8|] and the derivation of the unique non self-referential (/-canonical 

nn 

distribution 19] [20], are also based on the g-product. 

VI. CONCLUSION 

Starting from a fundamental nonlinear equation dy/dx = y q , we present the scaling 
property and the algebraic structure of its solution. Moreover, we prove that the algebra 
determined by its solutions is mathematically consistent with Tsallis entropy only as the 
corresponding unique information measure based on the following 2 mathematical reasons: 

(1) derivation of Tsallis entropy from the g-multinomial coefficient and g-Stirling's formula, 

(2) coincidence of the identity derived from the g-multinomial coefficient with the generalized 
Shannon additivity which is the most important axiom for Tsallis entropy. 

These mathematical discussions result in the self-consistency between the mathematics 
derived from the g-exponential and Tsallis entropy. 

Recently, we have presented the following fundamental results in Tsallis statistics: 

1. Axioms and the uniqueness theorem for the nonextensive entropy 1 

2. Law of error in Tsallis statistics 

3. g-Stirling's formula in Tsallis statistics jl^ 

4. g-multinomial coefficient in Tsallis statistics [U 

5. Central limit theorem in Tsallis statistics (numerical evidence only) j3| 

6. g-Pascal's triangle in Tsallis statistics ^| 
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7. The unique non self-referential g-canonical distribution in Tsallis statistics [19] [20] 

8. Scaling property characterized by the fundamental nonlinear differential equation [the 
present paper]. 

All of the above fundamental results are derived from the algebra of the g-product and 
mathematically consistent with each other. This means that the g-product is indispensable 
to the formalism in Tsallis statistics. More important point is that the g-product originates 
from the fundamental nonlinear differential equation (JTJ with scale-invariant solutions. 
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